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Abstract. We compare three different models of clumpy gas disk and show that the Circumnuclear Disk (CND) in the Galactic 
Center and a putative, geometrically thick, obscuring torus are best explained by a collisional model consisting of quasi-stable, 
self-gravitating clouds. Kinetic energy of clouds is gained by mass inflow and dissipated in cloud collisions. The collisions give 
rise to a viscosity in a spatially averaged gas dynamical picture, which connects them to angular momentum transport and mass 
inflow. It is found that CND and torus share the same gas physics in our description, where the mass of clouds is 20 - 50 M® 
and their density is close to the limit of disruption by tidal shear. We show that the difference between a transparent CND and an 
obscuring torus is the gas mass and the velocity dispersion of the clouds. A change in gas supply and the dissipation of kinetic 
energy can turn a torus into a CND-like structure and vice versa. Any massive torus will naturally lead to sufficiently high mass 
accretion rates to feed a luminous AGN. For a geometrically thick torus to obscure the view to the center even super-Eddington 
accretions rates with respect to the central black hole are required. 
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1. Introduction 

The central black hole of the Galaxy is surrounded by several 
hundred clouds of gas and dust forming a disk-like structure 
(Circumnuclear Disk CND, see e.g. Giisten et al. 1987) up to a 
radius of ^7 pc 1 . The CND has a hydrogen mass of a few times 
10 4 M Q which is distributed in clumps with an estimated area 
filling factor of <f>& ~ 0.1 and a volume filling factor of </>y ~ 
0.01. The clumps have densities of several 10 5 cm -3 , radii of 
~ 0.05 pc and gas temperatures > 100 K. A typical dusty clump 
has a visual extinction Ay > 30 m and Mh ^30 M Q . Further 
properties of the clouds in the central 2 pc are listed in Jackson 
et al. (1993). The vertical thickness (2 H) of the CND increases 
from ^0.5 pc at a radius of 2 pc to about 2 pc at the outer radius 
of 7 pc. The disk rotates with a velocity of ^100 km s _1 which 
corresponds to a Keplerian velocity around a central object of 
several 10 6 M Q and has a velocity dispersion of ~ 20 kms -1 . 
With an R/H of ~ 5 — 7 (R being the radius to the center), the 
CND can be qualified as a thin disk. At a distance of ~ 1.5 pc 
the molecular gas density drops by more than a magnitude, i.e. 
the central 3 pc of the Galaxy are void of molecular gas. In this 
sense the CND resembles more a ring or a torus. 
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1 We assume 8.5 kpc for the distance to the Galactic Center. 



In active galactic nuclei (AGN) the central luminous and 
compact object (e.g. the inner accretion disk around the black 
hole) is likely to be obscured by the torus, if the mean free 
path of clouds along the line of sight is greater than one. This 
qualitatively obvious statement was quantified by Nenkova et 
al. (2002) based on the comparison between models of cloudy 
dust emission and IR spectral energy distributions. They con- 
clude that the number of clouds along the line of sight should 
be 5-10, which corresponds to the mean free path of clouds 
being 5-10 times smaller than the radius of the torus. Krolik 
& Begelman (1988) were among the first to discuss the overall 
physical properties of the circumnuclear gas in active galactic 
centers. They already emphasized the importance of the cloudy 
structure of the circumnuclear material. 

Based on the unified model (Antonucci 1993) the spectro- 
scopic difference between Seyfert 1 and 2 nuclei is due to the 
line of sight passing through the torus (Sy 2) or not (Sy 1). 
From the relative numbers of Sy 1 and Sy 2, Schmitt et al. 
(2001) derive a half-opening angle of 48° for the obscuring 
torus. But Tran (2003) shows in a recent paper that about half 
the Sy 2 nuclei derived from a 12/im sample do not have de- 
tectable hidden broad line regions (HBLR) in spectropolarimet- 
ric observations. Because the obscuration determined from X- 
ray spectra is not systematically different for Sy 2 with HBLR 
and those without HBLR, it is not obvious that the obscuring 
material in non-HBLR Sy2 nuclei, which seem to be intrinsi- 
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cally fainter, is in form of a torus. It is therefore difficult to 
determine the half-opening angle of the torus and its relative 
thickness statistically. 

More evidence can be found by imaging extended NLR in 
Olll for example. For some sources the extended NLR take the 
shape of an ionization cone either due to a anisotropic radi- 
ation source or a collimating torus (for a review see Wilson 
1996, and Mulchaey, Wilson & Tsvetanov 1996 for a list of 
morphologically identified sources). The latter scenario is sup- 
ported by evidence for obscuration at the apex of some of the 
observed cones. If this torus is responsible for an ionization 
cone, it extends radially from a few pc to several tens of pc and 
its flaring should be a — (180° — ip)/2, where ip is the open- 
ing angle of the ionization cone. With a mean opening angle of 
ip = 80° ± 20°, the disk flaring angle is a ~ 50° ± 10° or 
H/R ~ 1, which means that these tori are geometrically thick. 

The CND in the Galactic Center in its current form is 
clearly not an obscuring torus, because the mean free path of 
the clouds along the radius is l c \ ~ 10 pc and thus R/l c \ < 1, 
where R = 7 pc is the extent of the CND, i.e. the line of 
sight through the torus. It is not a thick disk either, because 
H/R~ 0.2. 

In this article we investigate if there are common prop- 
erties between the CND in the Galactic Center and obscur- 
ing tori in Seyfert galaxies. In Sect. [2]we recall the theory of 
clumpy accretion disks. The equations needed for the compari- 
son between the different models and observations are given in 
Sect. [3] The models are than applied to the CND and a thick, 
obscuring torus (Sect. [4}- We discuss the results in Sect.|5]and 
give our conclusions in Sect. [6] 

2. The theory of clumpy gas disks 

2.1. Turbulent disks 

In a recent work (Vollmer & Beckert 2002, Paper I) we de- 
veloped an analytical model for clumpy accretion disks and in- 
cluded a simplified description of turbulence in the disk. In con- 
trast to classical accretion disk theory (see, e.g., Pringle 1981), 
we do not use the "thermostat" mechanism, which implies a di- 
rect coupling between the heat produced by viscous friction and 
the viscosity itself which is assumed to be proportional to the 
thermal sound speed. Thus, the viscosity, which is responsible 
for the gas heating, depends itself on the gas temperature. This 
leads to an equilibrium corresponding to a thermostat mecha- 
nism. Instead, we used an energy flux conservation, where the 
potential energy that is gained through mass accretion and dif- 
ferential rotation is cascaded by turbulence from large to small 
spatial scales. The inertial range of the turbulent cascade under 
consideration will not reach down to the scale of microscopic 
dissipation. Instead we regard the energy to be dissipated at 
an inner scale where the first cascade coming from the largest 
scales ends or breaks. 

In a first approximation this happens when the gas clouds 
become self-gravitating and dissipation then takes place inside 
individual clouds. Because the energy reservoir is the gravi- 
tational binding energy in the potential of the galaxy, a star 
cluster, or a central black hole, and the dissipation scale is 



set by self-gravity, we will call this scenario the fully gravi- 
tational model (FG model) in the following. Turbulence trans- 
fers the energy from the driving wavelength l^ liv to the dissi- 
pation wavelength Idiss, which corresponds to the length scale 
of a self-gravitating cloud. The two length scales are linked by 
the scaling parameter £ which measures the extent of the in- 
ertial range. For a Kolmogorov-like turbulent energy spectrum 

C = Cdriv/^diss) 3 - 

An additional 'free' parameter of the model is the Toomre- 
Q parameter: 



Q = 



VtUTb H 

ttGS 



> 1, 



(1) 



where Uturb is the turbulent velocity, SI the angular velocity, G 
the gravitational constant, and £ the gas surface density of the 
disk. For all models we assume Q to be constant throughout the 
disk. If we can approximate the total gas mass within a radius 
R by A / gas = hR 2 Yj Toomre-Q can be rewritten 



Q = 



VtuYb M tot 

^rot ^"eas 



(2) 



where M to t is the total enclosed mass and i> rot the rotation ve- 
locity. 

Like in a standard accretion disk scenario a viscosity al- 
lows redistribution of angular momentum and mass accretion 
together with a gain of kinetic energy from the potential well 
in which the disk resides. For the viscosity we use the following 
prescription: 



urb^drxv 



(3) 



A set of equations is obtained with 3 parameters: the Toomre 
parameter Q, the scaling parameter £, and the mass accretion 
rate within the disk M. These equations can be solved analyt- 
ically and the results were already used in Paper I to describe 
properties of our Galaxy's extended gas disk. Besides the pa- 
rameters above the solutions depend on the u rot , and R. Within 
the model it turned out that the driving wavelength in Eq.(|3} 
equals the disk height ^driv = H. 

In a second step we included the energy input due to super- 
nova (SN) explosions (Vollmer & Beckert 2003, Paper II). Here 
the energy flux provided by SNe is transfered by turbulence to 
smaller scales where it is again dissipated. The SN energy flux 
is assumed to be proportional to the local star formation rate p*, 
which, in turn, is taken to be proportional to the mean density 
and inversely proportional to the local free fall time of clouds. 
These clouds have sizes that are a factor 8 smaller than the driv- 
ing length scale. The factor of proportionality is the probability 
to find a self-gravitating cloud, i.e. the volume filling factor. 
Because the model works with vertically integrated quantities, 
like the surface density, the integration length in the vertical di- 
rection is assumed to be the turbulent driving length scale, i.e. 
the vertical length scale over which clouds can become self- 
gravitating: = p* /driv The SN energy per unit time EgN 
per area A A is therefore proportional to the local star formation 
rate £*: 



Esn 
~AA 



(4) 
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where the factor of proportionality £ is independent of the ra- 
dius in the disk. Its normalization with Galactic values yields 
£ = 4.6 10" 8 (pc/yr) 2 . 

In the FG model it is assumed that the energy flux trans- 
ported through the turbulent cascade is supplied by mass ac- 
cretion and the energy balance is 



pv 



y turb 



driv 



1 dn 

MVmt "r-r- 

2tt dR 



(5) 



In the case of SN driven turbulence the energy flux is deter- 
mined by the star formation rate 



pv 



y turb 



driv 



(6) 



In most gas disks on galactic scales the energy supply by su- 
pernovae dominates the gravitationally released energy for the 
same disk properties except for large Q values. The situation 
changes if the disks are so thin, that SN bubbles break out and 
most of their energy is lost from the disk. 

Furthermore, we take into account that the clouds are 
formed at the compressed edges during the interaction between 
SN remnants. The size of the clouds l c i is smaller than the tur- 
bulent driving wavelength and we use 6 = /driv/fci with 5 > 1. 

2.2. Collisional disks 

While the turbulent disk model describes the ISM as a con- 
tinuous medium where a turbulent cascade can develop, this 
changes when the self-gravitating clouds become stable. In 
such a medium an equilibrium disk can form, if there are frag- 
menting or partially elastic collisions (the clouds are magne- 
tized) between the clouds. These collisions give rise to a vis- 
cosity which depends on the velocity dispersion of the clouds 
(we will call it Uturb nonetheless) and the mean free path of 
cloud with respect to cloud collisions i co n. If the collisional 
timescale t co jj = Z C oii/^turb is larger or equal to the dynamical 
timescale, the resulting viscosity (Pringle 1981) can be written 
as 



v = t 1 v tUYh H , 



(7) 



where r = t co rjf2 and H is the disk height. Here it is already 
assumed that the disk is in hydrostatic equilibrium vertically 
H = vtmh I where the pressure is provided by the velocity 
dispersion. The energy dissipation rate due to collisions can be 
written as 



AE 
AA At 



2 2jt 'turb 
teal] 



2f^turb 
lean 



(8) 



The last equation on the right implies that the energy source 
is differential rotation and accretion as in the FG model. Here 
we introduced an efficiency parameter e which measures the 
fraction of the kinetic energy of clouds dissipated in a typical 
collision. Because most conclusions from this model are inde- 
pendent of the particular form of the energy equation ^ we 
set e = 2 /3 and so we are consistent with Eq. (0 and the use 
in Paper I. Since the driving wavelength in Eq. for the FG 
model is the disk scale height Zdriv = H, the FG model and 



the collisional model are formally equivalent when we identify 
t with £. But the interpretation of the viscosity in terms of a 
continuous turbulent medium (FG model) and the collisional 
viscosity here are completely different. 

3. A Comparison of characteristic properties of 
the three models 

In this section we recall relevant relations between observables 
within the FG and the SN models and derive observable prop- 
erties of the collisional disk. This will allow us later on to deter- 
mine the free model parameters and compare the models with 
observations. 

3. 1. The fully gravitational (FG) model 

Following Paper I this model has two dimensionless param- 
eters: Toomre-Q defined in Eq. (0 and £ in Eq. Q. From a 
solution consistent with observations it is possible to derive the 
mass accretion rate in the disk M from the FG model. Very im- 
portantly the turbulent velocity dispersion can be expressed as a 
function of the gravitational constant G, the scaling parameter 
£, the mass accretion rate M, and Q: 



Uturb = {\gCQM)~ 3 ■ 

Furthermore the cloud radius is given by 



(9) 



(10) 



and the volume filling factor, which is defined by 0y = p/pch 
where p is the mean density in the disk and p c \ is the density of 
the clouds, yields: 



1 



(ID 



For a disk with almost constant surface density the total gas 
mass is given by 



M gas = TtSi? 2 = 2-3G-3(3Q-3M3 Vrot JS 



(12) 



where the rotational velocity v TOt can be either directly ob- 
served or determined from the total mass distribution. In the 
case of an isothermal star cluster the surface density also has 
a cusp with S oc R^ 1 and the enclosed gas mass is a factor 2 
larger than in Eq. (I12> . 

3.2. The SN model 

In Paper II we distinguish two cases: (i) a self-gravitating gas 
disk in vertical direction (Q = 1) and (ii) Q > 1. 

For Q — 1 the velocity dispersion in the disk is given by 



w turb = 0.82GttM-<T-£tt 



(13) 



where £ is the normalization of the SN energy input in Eq. 
and S the ratio if actual cloud size to the driving length scale. 
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For Q > 1 the velocity dispersion in the disk and the driv- 
ing length scale are 



"turb = 0.87 X^GsMsQTS S~TE £15 ^ , 



(14) 
(15) 



where \ = H/R is the relative thickness of the disk. The vol- 
ume filling factor turns out to be 



0.79 x~~GsMsQ-ts 5~ 15 Vr js 



(16) 



It is possible to give a criterion for the validity of the model 
that is when the energy input through supernovae explosions 
Eq. l|6) dominates over the gravitationally released energy due 
to accretion (Eq.|5}. This can be turned into the following con- 
dition for the volume filling factor: 



-, 8 r 4 

y turb l driv 



(17) 



3.3. The collisional model 



For the collisional model we can get quite general relations 
for cloud sizes and a criterion for the Q parameter to be sat- 
isfied so that clouds are actually stable. As already mentioned 
in Sect. 12.21 the mean free path of a cloud is 



lr 



^turb 
trnU 



3urb 



tH 



(18) 



while the mean free path on the other hand is the inverse of the 
cloud number density n times the cloud cross section o = irr 2 ^. 



lr 



[ 4r c ipH> 



4r c i 
3 0v 



(19) 



In contrast to the turbulent models, the condition for clouds to 
be self-gravitating is that the free fall time equals the sound 
crossing time (in the case of the turbulent models it is the tur- 
bulent crossing time): 



ts = 



3tt 



32Gp c 



— ^sound — 



If! 



(20) 



where c s is the sound speed within the clouds. Inserting p = 
Vl 2 / (itGQ) from the definition of the Q-parameter (Paper I) 
leads to 



Q 



•r 1 r 2 

-'turb L s 



(21) 



Not surprisingly the size of clouds gets larger when the disk is 
further away from the gravitational instability (larger Q values) 
and the collision time is long. The critical density against tidal 
shear is p cr i t = 3f2 2 / (4irG) and thus p/ p cri t = 4/(3Q). On 
the other hand we have p/p c \ = <fiv- The condition for the 
clouds to resist disruption by tidal shear is p c \ > p cr it, which 
leads with the help of Eq. (I20> to an upper bound for the cloud 
size: 



r c \ < 



7T C s 

Tin 



(22) 



Table 1. Observed properties for the CND 



cloud radius 




0.05 pc 


cloud mass 




30 M 


sound speed within the clouds 


c s 


1.05 kms" 1 


turbulent velocity dispersion in the disk 


«turb 


20 kms -1 


disk radius 


R 


7 pc 


disk height 


H 


1 pc 


disk mass 


M sas 


210 4 M Q 


rotation velocity 


Viot 


115 kms" 1 



Inserting Eq. d2 II for the cloud size gives the following con- 
dition for the ratio of Toomre parameter to the dimensionless 
collision time r 



Q VS 



< 



Wturb 



(23) 



or using the definition of the Toomre-Q (Eq. ^ a lower bound 
on the collision time for stability of clouds 



r > 



(24) 



The upper bound on Q/t increases for thick disks (e.g. large 
velocity dispersions), while the limit on the collision time gets 
shorter for more massive disks. 

The limit on Q/t is especially interesting for geometri- 
cally thick and obscuring tori, because obscuration requires 
T = l co \i/H to be small, which turns into a severe constrain 
on the allowed Q-values. 

For the CND the clumpy structure is allowed to be transpar- 
ent, but with the known cloud sizes and temperatures (sound 
speed) it is possible to get a measure of r from the formally 
equivalent FG model expressions Eq. d9lll2> by replacing ( by 
r, and Eq. (ED 



R 2 



8 GM gas r cl 



(25) 



Obviously the collision time gets longer for smaller clouds, less 
massive disk, but surprisingly also for higher temperatures in 
the clouds. 

4. Finding the most appropriate model 

In this Section we will determine the model parameters of the 3 
different models for (i) the Circumnuclear Disk in the Galactic 
Center and (ii) for obscuring tori. This will enable us to discuss, 
which models are realistic for the two kinds of ring structures. 

4. 1. What it needs to make a Circumnuclear Disk 

The observed or from observations deduced properties of the 
CND are summarized in Tab. ^ The CND resides in a region 
where the central star cluster starts to dominate the total mass 
(e.g. Schodel et al. 2002). The total enclosed mass at radius R 
rises linearly with R and the rotation velocity is almost constant 
throughout the CND. Because the CND extends only over a 
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limited range in radius we use M gas = irHR 2 for the gas mass 
contrary to the discussion at the end of Sec. 13. II For all models 
the Toomre-Q parameter is 

„ VturbVrot-R 



Table 2. Assumed properties of an obscuring torus 



GA'h 



190 



(26) 



With this number fixed we can go to the individual models 



4.1 .1 . Massive clouds in the FG model 

From the cloud size and the disk height inserted in Eq. dlOt we 
get £ = 20. Together with the Q-value this leads to a volume 
filling factor V = {(Q)^ 1 ~ 2.6 10~ 4 . Independent of the 
volume filling factor we find from Eq. © a mass accretion rate 
of M = 10~ 3 M Q yr _1 . The mass of a single cloud follows 
M c i = (4/3) irr^p/cjjy = 260 M Q , but notice that this a rather 
strong function of the assumed outer radius of the CND, be- 
cause it enters quadratically in the mass density. The number 
of individual clouds in the CND is only 80 in FG model. 

4.1 .2. Why supernovae do not play any role in the 
Circumnuclear Disk 

In a first attempt we followed Paper II and assumed S = 5, but 
the resulting driving length scale was much to large. We could 
reduce that problem by choosing 5 = 2. The relative thickness 
of the disk is known \ = H/ R = 0.14, so that we derived from 
Eq. J14ft a mass accretion rate of M = 2.3 10 -3 M yr _1 . The 
expression for the volume filling factor (Eq. I16> . the driving 
length scale (Eq. 1 1 51 . and the cloud radius then give 0v = 
8.2 10~ 4 , Zdriv = 1-24 pc, and r cl = l dl - iv /S = 0.62 pc. In 
this scenario the cloud radius is derived from the model and 
is listed in Table for that reason. Consequently, the mass of 
a single gas cloud is M c \ = 1.6 10 5 M even larger than the 
whole mass of the CND. 

In a further approach we treated 5 as a free parameter, solv- 
ing Eq. ( 1151 and \14\ together with r c \ = ^driv/<5 leads to the 
incredible ratio 

5 = 0.86QC 5 — 4fe = 7.4 10 6 



(27) 



which states that ^ r ; v ps 370 kpc. 

In a last attempt we left £ the normalization of the SN 
energy input rate as a free parameter. This requires us to fix 
the driving length scale Zdriv, which up to now was always 
too large, to the disk scale height ^dnv = H = 1 pc. The 
cloud radius from Table Q is used as input here and is in- 
cluded in Table [3] for completeness only. This immediately 
leads to S = 20, an accretion rate of M = 6 10~ 3 M yr _1 , 
a volume filling factor of 0y = 10~ 6 , and a cloud mass of 
M c \ = 5 10 4 M . Apart from the much too large cloud mass 
one sees that this model also fails for the SN efficiency, because 
we find £ = 2.7 10 7 (pc/yr) 2 , which requires star formation 
and the resulting SN energy input into the disk to be about 600 
times more effective than elsewhere in the Galaxy. The large 
driving length scale or the large star formation/SN efficiency 
required shows, that energy supplied by SN under normal con- 
ditions is much too small to drive the dynamics of the CND. 



sound speed within the clouds 


c s 


1.2 kms" 1 


turbulent velocity dispersion in the disk 




72 kms- 1 


disk radius 


R 


10 pc 


disk height 


H 


6 pc 


rotation velocity 


"rot 


120 kms" 1 



4.1 .3. A transparent collisional CND 

The observed properties of the CND allow us to estimate the 
'transparency' t = l co \\/ H of the CND from Eq. d25i which 
results in t = 15. So on average there are about t j (2ir) w 2 
collisions for each cloud per orbit in this model. The expression 
for the velocity dispersion Eq. (|9j, which is equally valid for 
the collisional case with Q = t, gives an estimate of the mass 
accretion rate of M = 1.3 10~ 3 M Q yr _1 , which is almost the 
same as in the FG model, because t here and £ for in the FG 
case are not much different. The volume filling factor of clouds 
0v = (4r c ifl)/(3ru t urb) = 3.6 10 -3 is significantly larger 
than in the FG model. The reason is that Q enters in 0v of the 
FG model (Eq. II It . while it is the cloud size in the collisional 
case here. The result is a substantially smaller mass of typical 
clouds M c \ = 19 Mq. While the transparency derived from 
Eq. d25l > depends on global properties like the total gas mass in 
the CND, t can also be determined from a purely geometrical 
consideration 



n 



nav tulb 



20, 



(28) 



with the cloud number density n = N c \/V = (M gas /M c i) /V, 
the total number of clouds A^ cl = 1044, and V = 200 pc 3 the 
volume of the disk. This value is reasonably close to the one 
derived from Eq. d25i above. 

4.2. Building Geometrically and Optically Thick Tori 

For obscuring the BLRs in type 2 AGN and for collimating 
outflows and the ionizing flux in the extended NLR (ionization 
cones) we suggest properties for an obscuring torus as sum- 
marized in Tab. |5] The values would be appropriate for our 
galactic center becoming a Sy 2 nucleus with a compact dusty 
torus. The results are easily rescaled to objects like NGC 1068, 
where NIR-speckel observations (Wittkowski et al. 1998) show 
that the inner rim of the torus, if interpreted as such, is less 
than 2 pc from the AGN. The torus rotates somewhat faster and 
since there is an enhanced radiation from the active nucleus, 
the sound speed within the clouds is somewhat higher. 

As already mentioned in Sect.^the torus is obscuring if the 
mean free path of the clouds ^ co n is smaller than the line of sight 
through the torus. In a first approximation this corresponds to 
the condition / co n = H where we expect that in the midplane 
the line of sight is blocked by several clouds (Nenkova et al. 
2002) and the torus is also geometrically thick H ~ R. Within 
the framework of the FG and collisional models this leads im- 
mediately to the requirement £ = r = 1. 
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We can start the discussion with the statement that the FG 
model fails to describe a thick torus, because the requirement 
C = 1 implies that the cloud radius equals the disk height r c \ = 
H/( — H. Since this is unphysical, we can reject this model 
for obscuring tori. 

4.2.1 . Colliding clouds in a thick torus 

From the condition that the gas clouds have to resist tidal shear 
Eq. d23t reads Q < 54. As the torus is thicker for larger Q 
we assume Q to take its maximal allowed value. This means 
that typical clouds are actually at the shear limit. Together with 
t = 1 we find clouds of r c \ = 0.1 pc. From Eq. (|9j like in the 
CND case we find an accretion rate of M = 3.3 M Q yr _1 and 
from (Eq. (II 21 a gas mass in the torus of A/ gas = 4 10 5 Mq. 
The volume filling factor is determined in analogy to the CND 

~ 2 . The resulting mass 



case (Sec. 14.1.31 and gives cjty = 2.5 10" 
of a single gas cloud is here M c \ = 45 M Q . The total gas mass 
is larger than in the collisional CND case as expected, because 
more mass in clouds is needed for obscuration. The implied 
cloud sizes and cloud masses are surprisingly similar to those 
of the CND, which may indicate that the dominating physical 
processes in CND and torus are the same. The volume filling 
factor is of course larger as we would not get an obscuring torus 
otherwise. The higher density of clouds now lead to more cloud 
collisions and consequently the increased mass accretion rate, 
which feeds the AGN. 



Table 3. Derived parameters and properties of the CND and 
an obscuring torus (TORUS). The values in brackets for the 
SN model are for £ treated as a free parameter as described in 

Sec. EES 







CND 


TORUS 


units 


FG model 










Toomre parameter 


Q 


190 


— 




^driv/^diss 


C 


20 


1 




mass accretion rate 


M 


1 10~ 3 


— 


Moyr" 1 


volume filling factor 


4>v 


2.6 10~ 4 






cloud mass 


M c i 


260 




M 


SN model 
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4.2.2. Supernovae can drive tori at large mass 
accretion rates 

For the SN driven turbulent case of a torus model it is not ob- 
vious from the beginning what the appropriate Q-parameter is. 
The torus will only be geometrically thick for large turbulent 
velocities. This together with the normalization of the SN effi- 
ciency £ derived from integrated galactic values (Paper II), im- 
plies that a mass accretion rate (Ea.ll3> of M ~ 45 d 1 ^ 3 M Q /yr 
would be needed for a (Q=l)-torus. 

In the case of a torus with Q > 1 we can estimate, if the 
energy supplied by supernovae will dominate the heating rate 
due to viscosity and differential rotation. Inserting the volume 
filling factor from Eq. d!6i into Eq. J 1 71 leads to a limit for the 
product SQ 2 : 

SQ 2 < 44- ■ ( 29 ) 

With the values for the torus model from Table|2]it gives a limit 
5Q 2 < 8.6 C 3 - The comparison in the inequality d29t is made 
with the collisional case and for a torus it implies Q = t = 1. 
This results in a severe limit for both 5 and Q, because both 
parameters should be larger than 1. The range of allowed values 
for 6 and Q is rather small and we can assume (5 = 2 and 
Q = 2 for the moment. The model is now completely fixed 
and we find a huge mass accretion rate of M = 42 M yr _1 
through the torus, with cloud sizes of r c \ ~ 2 pc and cloud 
masses of 280 Mq. The torus in this model has a volume filling 
factor of clouds of 4>y ~ 6 10~ 2 and a total mass in gas of 



-A/g as ~ 1 10 7 Mq. The large cloud sizes and 5 = 2 implies 
that the driving length scale is l^ T i v = Q.7H almost as large as 
the torus scale height. Because Q is small in this model and the 
torus is required to be geometrically thick with H ~ R, the gas 
mass in the torus is comparable to the total enclosed mass in 
stars and central black hole. 

5. Prerequisites and Consequences 

In the previous section we have quantified what is needed for 
turbulent (fully self-gravitational, or supernovae driven) mod- 
els and a model consisting of quasi-stable clouds undergoing 
collisions, which lead to energy dissipation balanced by ac- 
cretion in the gravitational potential of the galactic center. A 
resume of the results for the Circumnuclear Disk and an ob- 
scuring torus is given in Tab. [3] 

For the Circumnuclear Disk (CND) in the Galactic Center 
the model of supernovae driven turbulence (SN model) can 
clearly be rejected, because the inferred cloud mass is larger 
than 10 4 M Q and all the mass in the CND would be only one 
cloud. Within the framework of the fully gravitational (FG) 
model, the volume filling factor is small (<f>y = 2.6 10 -4 ) and 
the cloud mass is about 8 times the observed value (M°^ s ~ 
30 M Q ). Even when allowing for the uncertainties in the ob- 
servational parameters like total mass and outer radius of the 
CND and the crudeness of the approximations in the model, it 
is nonetheless unlikely that this model is a realistic description 
of the CND. In contrast to the SN/FG-models the collisional 
model leads to disk and cloud properties that are close to the 
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observed values. This model thus gives the best approximation 
to the observations. It shares with the FG model the property 
that accretion is the energy source for keeping the velocity dis- 
persion, but it consists of stable, individual clouds, which are 
not part of a turbulent cascade. This conclusion is in qualitative 
agreement with those already drawn by Krolik & Begelman 
(1988). These authors derived the cloud properties on the ba- 
sis of an extrapolated, observed gas pressure. They then used 
a model for an equilibrium configuration where cloud mergers 
balance tidal shearing to show that the cloud area filling factor 
should be one. In this case the collisional time scale and the 
time scale for destruction by tidal shear are comparable. In our 
collisional model, the cloud properties are best described by the 
transparency r, which is not fixed to a singular value, and the 
fact that the clouds must resist tidal shear (Sec. l4.2t . 

5. 1. The efficiency of star formation and SN energy 
input in a collisional scenario 

In the following we compare the energy input due to SNe to the 
energy dissipation rate per unit area due to collisions within the 
collisional model described above. Supernovae can be impor- 
tant if their energy supplied to the ISM dominates the dissipa- 
tion rate due to collisions. This is nessesary but not sufficient 
to influence the dynamics of clouds, as we show below based 
on the ratio of cloud to SN remnant (SNR) surface area. 

For the first step of the argument the energy supply rate 
from Eq. @, modified by an efficiency 77, must be larger than 
the dissipation rate due to collisions. The star formation rate in 
a collisional scenario can be estimated following the discussion 
of Eq. (19) in Paper II by = 4>v^/ts- The free-fall time- 
scale is relevant for the size of clouds in the collisional case 
(Eq. I20i and r c \ taken from Eq. (121 > with Q at its maximum 
(Ea.l23> leads to the relation 



7 turb 

tH 



(30) 



From that we get a minimum efficiency r] e q for which SNe and 
accretion supply equal amounts of energy to be dissipated: 



'7c 



7T V. 



turb 



(31) 



For the collisional torus model described in Sec. 14.2.11 we find 
77 cq = 5.2 so that star formation and SNe are required to be 
more efficient than in galactic disks to become important. If 
the energy of SNe would be kept in the cloudy structure of the 
CND the equilibrium value would be ?7 cq .cND = 0.13, where 
we have used the 'observed' value for Q from Eq. i26\ . We con- 
clude that if star formation in the Galactic Center proceeds in 
the same way as in the Galactic disk, SN will provide a perhaps 
not negligible energy input in addition to dissipative collisions. 

Nevertheless, an equilibrium disk scenario where the dy- 
namics of clouds is driven by the energy input due to SNe is 
quite unlikely as we will show next. The SNR will expand into 
the inter-cloud medium (ICM). The ICM density and thermal 
velocity will limit the size the SNR and therefore the sphere 
of influence of the SN. To get an estimate for the ICM density 



we assume that it is in pressure balance with the clouds and 
fills the same region of space as the clouds do, which implies 
that the sound speed in the ICM equals the velocity dispersion 
Wturb- The ICM density is therefore pi CM = Pci c s/ U turb- Witr * 
the volume filling factor of the collisional model 

6 r 2 



it 2 Q 



(32) 



where p is the mean density in torus or CND. It turns out that 
the ICM density is 70% of the mean density for the CND and 
about 1% for the torus model; in both cases we find a parti- 
cle density of roughly 10 3 cm~ 3 , which is nice agreement with 
the result of Erickson et al. (1994). We can now estimate the 
maximum SNR size by using the power law fits provided by 
Thornton et al. (1998). To do that one has to realize that the 
SNR resolves when the shock speed of the SNR falls to the 
sound speed in the external medium which is the velocity dis- 
persion Wturb here. With v s cx we get from extrapolating 
the Thornton et al. result 



i? SN = 1.1 EH 1 V 42 



(Z/Z )-°- 1 pc. 



(33) 



In the CND most SNR are likely to break out of the relatively 
thin disk of molecular clouds and the efficiency will be largely 
reduced by this effect. In the case of the torus the fraction of 
SN energy impinging on clouds is simply the cross section of 
one cloud times the mean number of clouds within the SNR of 
radius Eq. J33I divided by the surface area of the SNR. This 
results in 77 = RsN't'v / (^ r c\) = 0y-RsN/(3r_ff) and the effi- 
ciency of SNe energy from star formation to drive cloud mo- 
tion is -q = 1.5 10~ 3 for the torus model. Thus we conclude 
that both CND and obscuring torus can be successfully de- 
scribed by the collisional model. The SN energy will almost 
exclusively heat the inter-cloud medium and they will have lit- 
tle effect on the dynamics of clouds in a scenario where clouds 
are stable and 'long-lived'. Krolik & Begelman (1988) used 
different arguments to estimate the energy input due to SNe 
and stellar winds, which is necessary to maintain a thick torus. 
They also found that the stirring of the circumnuclear mate- 
rial due to stellar driven processes (supernovae and/or stellar 
winds) is always weaker than the effects of cloud collisions. 

5.2. The transition from a CND to an obscuring 
torus and vice versa 

From the discussion up to now it seems that both, the 
Circumnuclear Disk in the Galactic Center and an obscuring 
torus, are well described by our collisional model. The den- 
sity of clouds within the model relative to the critical density 
against tidal shear p cr i t is p c \/ p cr it = 4/ (3Q</>v)- For the CND 
we find Pci/pcrit = 1-8 and the model for the obscuring torus 
is constructed in a way so that Q takes its critical value from 
Eq. d23t which automatically means that clouds are marginally 
stable against tidal shear. Both clump densities are close to the 
critical density with respect to tidal shear and we may spec- 
ulate that both structures are governed by the same physical 
processes. 

In this sense the CND would be an obscuring torus if it 
had more mass and it would be thick if it had a higher velocity 
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dispersion. We can quantify this with the definition of the (inte- 
grated) Toomre parameter Q (Eq.|2} to get = Q , where 
Mtot is the total enclosed mass and Af gas the enclosed gas mass 



at that radius. If clouds are indeed at the shear limit (Eci.l23> we 
can replace Q and get an expression for the transparency 



Aft, 



y/8 V rot M g 



(34) 



Because total mass and rotation velocity are given by the envi- 
ronment in the galactic center and the sound speed within the 
clouds is almost constant, we see that the torus becomes ob- 
scuring (that means r ~ 1) for large gas masses in the torus. 
Inserting Eq. J34t and the corresponding expression for Q in 
Eq. l|9) to derive the required accretion rate shows that this rate 
is M oc (iJ/i?)fturb-Afg as . Therefore large gas masses in a 
geometrically thick torus can easily provide the accretion rate 
needed to trigger and maintain nuclear activity. The Eq. J34l > 
for the transparency also shows that the torus around an AGN, 
which runs out of fuel and has less mass in in its surround- 
ing cloud distribution, will become transparent and eventually 
turn into a CND. Activity in a galactic center including a pos- 
sibly obscuring torus can also be restarted when additional gas 
masses in the form of molecular clouds arrive in the central 
region. 

With the critical gas mass for obscuration (r = 1) from 
Eq. J34i we can further see that the required mass accretion 
rate through the torus is 



Aft, 



M = —= I — ] c s 

V2 \RJ R 



(35) 



If the inner radius is set by the dust sublimation radius (see be- 
low) and the temperature in clouds is determined by radiative 
equilibrium with the central source, the sound speed is constant 
at the inner radius independent of the accretion rate or lumi- 
nosity. With (R = i?i n oc M 1 / 2 ) we find that the minimum 
accretion rate for obscuration is 



Al 



obs 



3.1- 



M© 
yr 



V- 



1/3 2/3 



(I) 4 ' 3 


c s 


2/3 




1.5km/s 





(36) 



where rj-\ is the radiation efficiency at the black hole rel- 
ative to 10% efficiency and Mr the enclosed mass in units 
lO 7 A/ . In addition it is assumed that all the mass transported 
through the torus actually reaches the black hole. The required 
mass accretion rates are highly super-Eddington MEdd = 
0.22M B H,7Af Q /yr for the central black hole with 10% effi- 
ciency for transforming rest mass energy into radiation. The 
solution is that most of the mass accreted through the torus is 
not effectively turn into radiation at the black hole. It is either 
blown away in an outflow or jet before reaching the black hole, 
and thus constitutes a cooling process for the remaining matter, 
or the accretion process is indeed radiatively inefficient at these 
super-Eddington rates. 

Since these collisional disks are stable and long-lived, 
a possible formation mechanism of a thick, obscuring torus 
might be the infall of large gas clumps onto a pre-existing cir- 
cumnuclear ring as investigated by Vollmer & Duschl (2002). 
Infalling gas clouds with different angular momentum will also 



lead to increased velocity dispersion in the merged system as 
required for a torus. 

5.3. The physics of clouds and the inner boundary 
of a torus 

In Vollmer & Duschl (2001a, b) we investigated the physics of 
gas clouds in a turbulent or collisional, clumpy disk. It turned 
out that under realistic conditions, where the gas is heated by 
the UV radiation of the central Hel star cluster, the gas clouds 
are stable and self-gravitating. Their radius is mainly deter- 
mined by the location of the ionization front due to the same 
UV radiation. Under realistic condition this results in the ob- 
served cloud size of ~0.1 pc. The present study suggests that 
the gas and dust physics of thick, obscuring tori is of similar 
nature. 

In this picture two selection effects operate on the initial 
cloud mass distribution: (i) only clumps that are dense enough 
to resist tidal shear can survive and (ii) clumps that are more 
massive than the Jeans mass collapse. In addition, the central 
UV radiation field determines the size of the clumps, i.e. their 
mass for a given central density. For the CND it turned out 
that at a critical distance of ~1.5 pc from the Galactic Center 
the clouds, which are dense enough to resist tidal shear, will 
collapse. This is in agreement with the observed inner edge 
of the CND. If such a mechanism is also at work in a thick, 
obscuring torus around an AGN, the inner edge of this torus 
would be determined in the same way, if the dust sublimation 
radius and the cloud evaporation radius are smaller than the 
critical radius for collapse and shear (C&S-radius) described 
here. 

Dust sublimation due to the radiation of a central source de- 
termines the size of the inner dust free cavity, which will set the 
source size for NIR continuum observations (e.g. Wittkowski et 
al. (1998) for NGC 1068). This is also the inner radius where 
cool clouds are shielded by dust opacity. A possible description 
of the sublimation radius by Dopita et al. (1998) 



R suh = 0.84pc 



sub 



1200K 



j 1/2 -1/2 
-^44 "-5 



(37) 



shows that f? sl ,b is a strong function of the sublimation tem- 
perature T su b, which depends on the dust chemistry and grain 
sizes with a_5 the typical grain radius in units of 10~ 5 cm. The 
determination of the sublimation radius assumes radiative equi- 
librium in the dust sublimation region and is due to the R~ 2 - 
dilution of the flux from the central luminosity L44 in units of 
10 44 erg/s. 

If the dust sublimation radius is larger than the C&S-radius 
described above, clouds with masses larger than ~10 M Q can 
only survive up to this radius. Inside, only the stripped cores 
of such clouds can survive (the "light clouds" in Vollmer & 
Duschl 2001a), because the location of the ionization front is 
determined by self-shielding and not dust absorption. 

If the cloud evaporation radius is larger than the C&S- 
radius, clouds evaporate before they can reach the C&S-radius 
or dust sublimation radius. Following Rrolik & Begelman 
(1988) clouds will evaporate when the evaporation time scale 
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equals the accretion time scale at a given distance from the ion- 
izing source, which leads to at a typical evaporation radius of 

-Revap 0.79 M gaS) 5 Z44 

= — TT? 1 (JO) 

PC M torus R out N a ,24 T 5 

for a cloud with a gas column density in units of 10 24 cm~ 2 , 
a temperature at the sonic point of the evaporating flow 
T5 of 10 5 K, a central ionizing luminosity L44 in units of 
10 44 ergs -1 , the mass accretion rate Mtorus measured in 
M Q /yr, the total gas mass in clouds M gaS) g relative to 10 5 M Q , 
and the outer radius of the torus i? out in parsec. Here we as- 
sumed E oc J? -1 , which is expected for a constant rotation 
curve (the torus residing in the potential of an isothermal star 
cluster). Inserting the values for the collisional torus model 
from Table[3]one finds i? C vap = 0.1 pc. The evaporation radius 
is sensitive to the ionizing luminosity and Nc\ in an individual 
cloud, so that Eq. J38i can be used to set a lower limit to the 
column density at a given radius. Eq. J38I seriously underesti- 
mates the evaporation radius if it gives i? C va P > -Rout/2. We 
can use this in the collisional model discussed here for a cen- 
tral source radiating at the Eddington limit of a 3 10 6 M© black 
hole and find that only clouds with Nq\ > 2 10 22 cm -2 will 
survive when they are exposed to the central radiation source. 
The clouds in the model with rci = 0.1 pc and 45 M are 
safely above that limit. At an inner edge of about 2 pc the 
clouds have a column density of > 10 24 cm~ 2 (assuming <fiv a 
constant) and all clouds in the model are stable against evapo- 
ration. 

6. Conclusions 

We explored three distinct models of clumpy gas disks for the 
Circumnuclear Disk (CND) in the Galactic Center and for a ge- 
ometrically thick, obscuring torus, as suggested in unified mod- 
els of AGN. Three kinds of models were tested: (i) a turbulent, 
fully gravitational model, where the energy, which is gained by 
mass inflow and differential rotation, is cascaded by turbulence 
to smaller scales. Alternatively (ii) the turbulence may be main- 
tained by the energy input due to supernova explosions, and we 
investigated (iii) a collisional model, where the gas viscosity is 
due to partly dissipative cloud collisions. In the cases (i) and 
(ii) we identify the energy sink to be the self-gravitation of the 
gas clouds, which separates them from the turbulent cascade. 

We find that the fully gravitational model cannot be entirely 
excluded for the case of a CND. The cloud masses for the CND 
model are larger than indicated by observations, but the masses 
in the model depend strongly on the assumed outer radius for 
an equilibrium disk structure. A turbulent cascade may also 
exist in a massive torus at exceedingly large mass accretion 
rates, where supernovae can drive the turbulence. Nonetheless 
the CND and a thick, obscuring torus are best described by the 
collisional disk model. We suggest that both the CND and a 
thick, obscuring torus share the same gas physics: they consist 
of stable, self-gravitating clouds that are heated by the central 
UV radiation and the density of clouds is in both cases close 
to the critical density with respect to tidal shear. The mass of 
single clouds is a few times 10 M Q in the CND and in a torus. 



In this picture a thick, obscuring torus is a more massive 
CND-like structure, that is with more gas clumps, and in addi- 
tion a higher velocity dispersion. We show that a torus provides 
automatically the high mass accretion rate needed to trigger and 
maintain nuclear activity and that a geometrically thick torus is 
obscuring the view to the center only at super-Eddington rates. 

The model for clumpy collisional disks suggests that the 
appearance as a CND or an obscuring torus only depends on 
the mass supply from larger radii in the galaxy, possibly in the 
form of pieces of giant molecular clouds. 

Acknowledgements. Part of this work was support by the DFG 
through grant SFB439 (WJD). 

References 

Antonucci, R. 1993, ARA&A, 31, 473 

Dopita, M. A., Heisler, C, Lumsden, S., Bailey, J. 1998, ApJ, 498, 
570 

Erickson, E. F, Colgan, S. W. J., Simpson, J. P., Rubin, R. H., Hass, 
M. R. 1994, Far-Infrared Line Observations of Sgr A West. In The 
Nuclei of Normal Galaxies, ed. R. Genzel and A. I. Harris, NATO 
Advanced Science Institutes Series C, Vol. 445, 249 
Glisten, R., Genzel, R., Wright, M. C. H. et al. 1987, Apl 318, 124 
Jackson, J. M, Geis, N., Genzel, R. et al. 1993, ApJ, 402, 173 
Krolik, J. H. & Begelman, M. C. 1988, ApJ, 329, 702 
Mulchaey, J. S., Wilson, A. S., & Tsvetanov, Z. 1996, ApJ, 467, 197 
Nenkova, M., Ivezic, Z., & Elitzur, M. 2002, ApJ, 570, L9 
Pringle, J. E. 1981, ARA&A, 19, 137 

Schmitt, H. R., Antonucci, R. R. J., Ulvestad, J. S., Kinney, A. L., 

Clarke, C. J., Pringle, J. E. 2001, ApJ, 555, 663 
Schodel, R., Ott, T, Genzel, R., et al. 2002, Nature, 419, 694 
Tran, H. D. 2003, ApJ, 554, L19 

Thornton, K., Gaudlitz, M., Janka, H.-Th., Steinmetz, M. 1998, ApJ 
500, 95 

Vollmer, B. & Duschl, W. J. 2001, A&A, 367, 72 
Vollmer, B. & Duschl, W. J. 2001, A&A, 377, 1016 
Vollmer, B. & Duschl, W. J. 2002, A&A, 388, 128 
Vollmer, B. & Beckert, T. 2002, A&A, 382, 872, (Paper I) 
Vollmer, B. & Beckert, T. 2003, A&A, 404, 21, (Paper II) 
Wilson, A. S. 1996, VA, 40, 63 

Wittkowski, M., Balega, Y., Beckert, T, Duschl, W. J., Hofmann, K.- 
H, Weigelt, G. 1998, A&A, 329, L45-L48 



